typeset using JPSJ.sty <ver.l.Ob> 



Kondo Effect in Single Quantum Dot Systems 
Study with Numerical Renormalization Group Method 



Wataru IzumidaQ, Osamu Sakai and Yukihiro S 



HIMIZU 



Department of Physics, Tohoku University, Sendai 980-8578 
^Department of Applied Physics, Tohoku University, Sendai 980-8579 

(Received March 16, 1998) 



The tunneling conductance is calculated as a function of the gate voltage in wide temperature 
range for the single quantum dot systems with Coulomb interaction. We assume that two 
orbitals are active for the tunneling process. We show that the Kondo temperature for each 
orbital channel can be largely different. The tunneling through the Kondo resonance almost 
fully develops in the region T <, O.lTjJ ~ 0.2T^, where is the lowest Kondo temperature 
when the gate voltage is varied. At high temperatures the conductance changes to the usual 
Coulomb oscillations type. In the intermediate temperature region, the degree of the coherency 
of each orbital channel is different, so strange behaviors of the conductance can appear. For 
example, the conductance once increases and then decreases with temperature decreasing when 
it is suppressed at T = by the interference cancellation between different channels. The 
interaction effects in the quantum dot systems lead the sensitivities of the conductance to the 
temperature and to the gate voltage. 
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§1. Introduction 

Coulomb oscillations behavior in quantum dot systems is well known as one of the typical phe- 
nomena due to the Coulomb repulsion between electrons in the dot. But it is one aspect of the 
electron-electron interaction effects. Local spin moment is induced by the Coulomb repulsion in 
the dot sometime. It will fluctuate by exchanging spin moment with leads and thus cause the 
strong inelastic scattering. But the local spin disappears in the lowest temperature due to the 
spin singlet formation between leads and the dot. These effects have been investigated as the 
Kondo effect for the magnetic impurity systems in many years ,0) and for the quantum dot systems 
recently.&i&iilii0§) 
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The origins of the non-vanishing spin moment in the dot are not unique. Of course unpaired 
electron in the odd number electrons in the dot brings the spin moment. Furthermore, degenerate 
orbitals, or strong Hund's rule coupling can induce the non-vanishing spin moment even in the even 
electron number cases.0^ This will cause varieties in the behavior of the spin state when the electron 
number is changed by the gate voltage. In the low temperature range that the Kondo singlet state 
is formed, the inelastic scattering due to the spin fluctuation is suppressed. Therefore the coherent 
process will dominate the tunneling. This process is considered as the resonant tunneling through 
the Kondo peak at the Fermi level in some sense. The Kondo peak sensitively diminishes as the 
temperature increases and disappears at high temperatures because the characteristic temperature 
for the many body singlet formation, the Kondo temperature, is usually very low. In general, the 
tunneling conductance will show various types of the interference effect depending on the geometry 
and the connectivity of the dot systems. The coherency of the system is disturbed by the inelastic 
scattering due to the spin moment revived at finite temperatures. However, calculation of the 
conductance in such the intermediate temperature has not been done theoretically because the 
reliable method to treat the Kondo effect has not been established. 

In this paper, we consider the system shown in Fig. 1, a quantum dot and connected two leads. 
We calculate the tunneling conductance in wide temperature range from the high temperature 
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Fig. 1. The configuration of the system. The z-axis is perpendicular to the xy plane. 

region in which the usual Coulomb oscillations are observed, to the lower temperature region in 
which the Kondo singlet state is formed. Conductance is calculated as the function of the gate 
voltage. We assume that two orbitals of the dot are active for the tunneling process. The following 
two cases are considered for the tunneling Hamiltonian (TH) between the conduction channels in 
the leads and the orbitals in the dot: (THi) There is one conduction channel in each lead, and it 
hybridizes with both of two orbitals in the dot. This means that the tunneling processes via the 
two orbitals interfere each other. (THii) There are two conduction channels in each lead, and each 
of the channels hybridizes separately with one of two orbitals in the dot. In this case the tunneling 
processes of the two channels are independent and do not interfere. (But the electronic state of the 
dot itself is affected by the interaction between electrons in the different channels.) In such two 
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situations, we study the following three cases for the orbital state in the dot; (1) orbital energies 
in the dot are not degenerate, (2) orbital energies in the dot are degenerate, (3) electrons on the 
orbitals have strong Hund's rule coupling. The six combinations of the situation, from THi-1 to 
THii-3, are studied systematically. 

There have been several methods to calculate the conductance in finite temperatures; perturba- 
tion theory, non-crossing approximation, quantum Monte Cairo method, numerical renormalization 
group (NRG) method.!' IBi'BS' The NRG is a reliable method in the temperature region T <, Tk, 
where Tk is the Kondo temperature.0'lll'0© This method directly calculates the tunneling ma- 
trices between wave functions of many electron states. Therefore we can obtain the tunneling 
probabilities in a general way for the coherent process in the low temperatures and the incoherent 
process in the high temperatures. The numerical calculations are performed by using NRG method 
through this paper. 

We found that the tunneling through the Kondo resonance fully develops in the temperature 
region, T ^ O.lTj^ ~ 0.2Tj^, where Tj^ is the lowest Kondo temperature of the dot system when 
the gate voltage is varied. The interference effect of the tunneling process also becomes important 
in this region. In high temperature region T ^ T^, the Coulomb oscillations type conductance is 
observed commonly. In the intermediate temperature range among these, various types of conduc- 
tance behaviors appear, especially in THi cases, because the coherency of each orbital channel is 
sensitively changed by the temperature and the gate voltage. 

In §^, the cases classified to the THi are studied, and the cases classified to the THii are studied 
in §^. In summary and discussion of this study are given. 

§2. One Conduction Channel in Each Lead 

To study the system shown in Fig. 1, following things are contained in the model: i) Electronic 
states in the two leads are written as the conduction band states, ii) There are two orbitals in 
the dot. iii) Electron-electron interactions in the dot are considered. And, iv) the electrons have 
tunneling matrices between the leads and the dot. For simplicity, we consider the situation that 
there are only two orbitals in the dot instead of considering many orbitals. We assume that the 
system has the mirror symmetry with respect to the yz and zx planes. (The z-axis is perpendicular 
to the xy plane.) 

In this section, we discuss the situation classified to the THi that there are only one conduction 
channel in each lead, and tunneling processes via the different orbital channels interfere each other. 
We consider the one-dimensional degree of freedom in each lead. This situation would appear in 
the case that leads are very narrow, or the joint part between the leads and the dot are very narrow. 
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2.1 Model Hamiltonian 

We consider the two orbitals in the dot, which are even and odd, respectively, under the yz 
mirror, and are even under the zx mirror. These orbitals in the dot are named as the even orbital 
and the odd orbital, respectively. We consider the one-dimensional degree of freedom with even 
symmetry under the zx mirror in each lead. We use the following Anderson Hamiltonian, 

H = Hi + Hd + Hi_d, (2.1) 

Hl = Yl Ekclk^CLka + Yl ^k4ika(^Rka, (2-2) 
k(T k<T 

Ha = (ea - XI + (^d + *) X! ^o<^ 

+ <^H ('^)cri<T2 ■ ('^)<T3<T4<^eiTl'^eiT2'^0(T3'^0(T4) (2-3) 

IT1<T2IT3IT4 

-f^l-d = -J^Y^ K4<t (cLfc(T + CRfea) + X! ^o'^oa (cLfca " CRfca) + h.C. (2.4) 
^ ka ^ ka 

The terms H\ and give the electron orbitals in the two leads and in the dot, respectively. 
The term i?i_d gives the electron tunneling between the two leads and the dot. The annihilation 
operator of the left (right) lead state is denoted by C]_,ka {cRka)^ and that of the even (odd) orbital 
in the dot is denoted by d^a {doa)- The quantity is the energy of the lead state. The quantity 
£d corresponds to dot's potential, and is able to change by applying gate voltage. The energy level 
of the even (odd) orbital is given by £d — * (^d + *)• The energy separation between the even and 
odd orbitals is defined as Ae^ = 2t. The Coulomb interaction constant is given by U for both the 
intra- and inter-orbital terms. We consider the exchange interaction between electrons on the even 
and odd orbitals with strength Jh- The operator cr is the Pauli matrix. The quantity (K)) is 
the hybridization matrix for the even (odd) orbital. (We neglect the k dependence of the tunneling 
matrices.) 

We introduce the even and the odd combinations of the lead orbitals, s^^ = {c]_,ka + CRka)/V^ 
and ajfco- = (cLfeo- ~ c^ka)/V^, respectively. The terms Hi and -ffi_d are rewritten as. 



Hl = Y ^ksl^Ska + Y ^k^ka'^ka, (2-5) 
ka ka 

Hi-d = E {Vedlska + h.c.) + Y {Vodlaka + h.c.) . (2.6) 

ka ka 

The hybridization strength for the even (odd) orbital channel is denoted as Ae = TrV^p{ep) (Aq = 

7rV^^p(£F)), where /9(£f) is the density of states of the lead states. Hereafter we denote the even 
and the odd orbitals by p as p = e and p = o, respectively. 
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2.2 Linear response conductance 

In this paper we restrict ourselves to the linear response conductance for the applied bias voltage 
between the left and the right leads. We generally define an electric current from the left lead L to 
the right lead R as follows; 

where — e is the charge of an electron, (A^l) is the expectation value of the time differentiation of 
the electron number operator in the lead L. The quantity {N^) is that of the lead R. We obtain 
the following expression for the linear response conductance formulajll^ 

= — lim 2.8 

h oj-*0 uj 



with 

4 Z 



m) 



x6 {lo - (En - Em)) , (2.9) 
where Z = J2n is the partition function, /3 is the inverse of the temperature, /? = 1/T. (See 



Appendix of the ref. 10 for the derivation of eqs. (2.8) and (2.9).) We note that we do not have done 
restrictive assumptions for the electric current such as the 'coherent tunneling' or the 'sequential 
tunneling' in the derivation of (2.8). The conductance for the case of THi, G*-^-*, is calculated by 
using eqs. (2.8) and (2.9) with A'l = Y.ka A.ka'^^ka and TYr = Y.ka ^^-Rku^^ka ■ We note that the 
tunneling processes via the even and the odd orbitals interfere each other in the tunneling between 



leads. Following the method described in ref. 10, the current spectrum, P (w), is calculated by 
NRG method. 

Here we briefly consider the tunneling at absolute zero temperature. The system would be the 
local Fermi liquid state at very low temperatures. In this case there are no spin scattering processes. 
In such the coherent tunneling process at zero temperature, the conductance for the case of THi, 
Gp\ is written as follows;^!!^ 

4^) = ^sin2||((ne)-(no))|. (2.10) 

This expression is the extension of the relation G-p = (2e^//i) sin^ {(7r/2) (n)} for the single orbital 
case.Bil' The quantity {n^) — (no) appears through the interference of the tunneling processes via 
the even and odd orbitals. In derivation of eq. (2.10), we have applied Friedel's sum rule to the 
even and the odd orbital channels separately.lll'0'0) 
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2.3 Numerical results with NRG method 
2.3.1 N on- degenerate orbital case 

In this subsection we present the numerical results for the non-degenerate orbital case. The 
parameters are chosen to satisfy the relations Ac, Aq ^ Aea: Ac = O.OOSvr, Aq = 0.0027r, Aej = 
0.1, [/ = 0.1 and Jh = 0. (In this paper we choose the band width D as an energy unit. The density 
of states for the conduction bands in each orbital channel are assumed to be constant within the 
region —1 < uj/D < 1.) 

First we present the electron occupation numbers at absolute zero temperature, (up) , as a function 
of the dot's potential, ^d- The conductance Gp^ at T = 0, which is calculated from (2.10), is also 

(2) 

shown in Fig. 2. (The conductance Gp at T = for the case of THii is also shown in the same 
figure. This quantity will be discussed in §p.) As the potential Ed is gradually decreased by applying 
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Fig. 2. The expectation value of the electron occupation numbers (up) for the even {p — e) and the odd (p = o) 
orbitals, and the conductance Gp"' for the THi case {n — 1) and the THii case (n = 2), as a function of the potential 
£d at T = for the non-degenerate orbital case. 



gate voltage, the first electron occupies the even orbital at about ~ ^ = 0.05 ((ric) ~ 0.5), and 
then the second electron occupies at about ^ t — U = —0.05 ((?ic) ~ 1.5). The third and 
the forth electrons occupy the odd orbital at about ~ — ^ — 2[/ = —0.25 ((no) ~ 0.5) and at 
8^ = —t — 3U = —0.35 ((no) ~ 1.5), respectively. (The potential values = t,t — U, —t — 2U and 
—t — 3U are just the cross points of energies of the states with different electron number in the 
dot, if one neglects the hybridizations between leads and dot.) The conductance Gp ^ has the finite 
value ~ 2e'^ /h (unitarity limit value) in the region oi t — U ^ Ed ^ t in which the even orbital is 
almost half-filled, and in the region of —t — 3U ^ £d ^ —t — 2U in which the odd orbital is almost 
half-filled. Both cases have odd electron numbers. Therefore, there will appear spin moment when 
the temperature rises near the Kondo temperature. This will cause inelastic scattering process and 
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will strikingly change the conductance from that of T = 0. 

Next, we present the temperature dependence of the conductance in Fig. 3. It is calculated from 
the current spectrum P (cj) by using eqs. (2.8) and (2.9). Let us see Fig. 3 from high temperature 




0.10 0.00 -0.10 -0.20 -0.30 -0.40 



Fig. 3. The temperature dependence of the conductance as a function of the parameter for the THi-1 case. (For 
the non-degenerate orbital case.) The regions given by broken hues are considered to be the magnetic state and 
sohd hnes region are the non-magnetic state (spin singlet state). This classification is done by comparing the 
temperature and the Kondo temperature estimated from the magnetic excitation spectra. The suffixes for the 
temperature are chosen to be common throughout § ^.3.1 and Fig. 15 in 



side. There are four peaks at the potentials ~ 0.05, —0.05, —0.25 and —0.35 in high temperature 
region T ^ 10"'^. These structures correspond to the usual Coulomb oscillations, because the each 
potential coincides with the point at which the electron number in the dot changes. The thermally 
broadened peak structures in Tio = 1.1 x 10^^ case change to more sharp peak structures in 
Tg = 3.8 X 10^^ case. When the temperature gradually decreases below about T ~ 10"^, the 
conductance grows up in the regions t — U ^ £d ^ t and —t — 3U ^ e^^ —t — 2U. The pairs of 
the peaks at high temperatures in the each region gradually change to the one peak structures. In 
the lowest temperature limit T ^ 10^^, the conductance G^^^ coincides to Gp^ in Fig. 2 except 
the numerical errors brought by the calculation of P" {uj)/uj with NRG method. (From comparison 
of the conductance g[!^ in Fig. 2 and G^^) at Ti = 5.8 x lO"'^ in Fig. 3, G^^) has about 15% 
smaller magnitude at = —0.30. This accuracy of the conductance through the calculation of 
P" {uj)/uj is not so good in a quantitative sense, but seems to be sufficient to extract the qualitative 
temperature dependence caused by the change of the electronic state.) As shown later, the Kondo 
temperature Tk strongly depends on ej- Hereafter, we denote as the lowest Kondo temperature 
when the potential Ed is changed. It is estimated to be about = 2.7 x 10^^ for = —0.30 case, 
which is near T4 = 1.6 x 10^^ in Fig. 3. The valley of the conductance at = —0.30 becomes 
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very shallow in the cases T ^ T3 ~ 0.2T^. The conductance changes from the usual Coulomb 
oscillations type to the Kondo resonance type at about T ~ 0.2T^. In the non-degenerate orbital 
case of this subsection, the interference between tunneling processes via the even and odd orbital 
channels seems to be not so operative. We expect that the temperature dependence of the growing 
up of the pair peaks to the fiat one peak is almost identical to that of the single orbital case. 

Next we present the excitation spectra in detail. Hereafter in this subsection, we mainly discuss 
only the region < —3U/2 = —0.15. In this case the Kondo effect of the odd channel electron 
is important. We have calculated the single particle excitation spectra for the even and the odd 
orbitals, Pp(a;) (p = e, o), which are given as follows, 



■l3Er, 



-/3E., 



X (^\{n\dl\m)\^ 6 {uj - {En - E^)) 
+ \{n\dp\m)\^6{u + {En-Em))). 



(2.11) 



The single particle excitation spectra at several temperatures for the case of = —0.30 are shown 
in Fig. 4. In Pe{io) we have no fine structures near the Fermi energy because the even orbital is 
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Fig. 4. The single particle excitation spectra at several temperatures for the potential ea = —0.30 ((rio) ~ 1.0, in 
which the odd channel is the Kondo regime). The results are for the non-degenerate orbital case. The spectra are 
plotted as a function of the logarithm of the energy except the inset plotted by the energy. 



almost full-filled. On the other hand the odd orbital is almost half-filled, (no) ~ 1.0, Po('^) has fine 
structures at low temperatures. The Kondo peak on the Fermi energy gradually grows up in the 
energy region ^ 10^"^ as temperature decreases from T 10^^. Growing up of the Kondo peak 
completes at about T <, 10~^. Reflecting these facts, the conductance near ~ —0.30 increases as 
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temperature decreases. The odd channel seems to be in the Kondo regime. In Fig. 5 the potential is 
raised to = —0.26 so that the odd channel is in the mixed valence regime ((rio) ~ 0.8). In Poi^^), 
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Fig. 5. The single particle excitation spectra at several temperatures for the potential ed ~ —0.26 ((no) ~ 0.8). 
The results are for the non-degenerate orbital case. In this case the odd channel is in the mixed valence regime. 



we can see an excitation peak near lu ~ —0.01 in the case of Tiq = 1.1 x 10^^. This corresponds 
to the excitation from the Nq = 1 state to the A'o = state. As the temperature decreases 
gradually, the peak position shifts to the Fermi energy, and the peak height increases. It grows 
up completely in the temperature region T ^ 10~^, and has structures seem to be a composite of 
the two components: the electron excitation of A'o = 1 — > A'o = and the Kondo peak. As seen 
in Fig. 4, these two components are separated with each other when the system is in the Kondo 
regime. Any way, the intensity of Po{'^) on the Fermi energy increases as temperature decreases, 
and causes the growth of the conductance in the potential region —t — 3U ^ Ed ^ —t — 2U. 

As seen from comparison of Fig. 4 and Fig. 5, the Kondo temperature is different in each 
potential case. To see this point clearly we calculate the magnetic excitation spectrum at zero 
temperature, Xm\^)^ 

x;M = EEK^I(^c,. + 5o,.)|Gr)|2 

n Gr 

Xd {to - (En - Ecr)) , (2.12) 

where Sp^z = {dp^dp^^ — dp^dp^i)/2 is the spin operator on the p-th. orbital and Gr denotes the 
ground state of the system. The peak position of the magnetic excitation reflects the characteristic 
energy of the spin fluctuation. The Kondo temperature Tk is defined as the energy of the peak 
position.Sil* It is shown in Fig. 6 as a function of the potential £d, and has the sharp local 
minima at = —0.30 (where (no) ~ 1-0) and = 0.00 {{ric) ~ 1.0). These minima characterize 
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0.10 0.00 -0.10 -0.20 -0.30 -0.40 
Fig. 6. The Kondo temperature as a function of the potential Sd for the non-degenerate orbital case. 

the temperature at which the conductance behavior changes from the Coulomb oscillations type to 
the Kondo resonance type. 

We point out that the Kondo temperature is also very sensitive to the change of the hybridization 
A as seen from comparison of the TK(ed = -0.30) = 2.7 x 10"^ (Aq = 0.0027r) and rK(ed = 
0.00) = 2.7 X 10"^ (Ac = O.OOStt), or from the analytical expression for the Kondo temperature. 
(Tk = \/UA/2 exp [—{ttU)/{8A) + A/{2U)] for the single orbital impurity Anderson model in the 
electron-hole symmetric cas ). This means that the Kondo temperature of the quantum dot 
systems is sensitive to the effective distance between leads and dot. 
2.3.2 Degenerate orbital case 

Next, we consider the situation that the two orbitals have degenerate energy levels. Even if the 
orbital energies are not degenerate in the strict sense, they should be considered as to be degenerate 
when the energy difference A^d is less than Tk. (See in Appendix.) The effect of spin fluctuation 
will appear even when electron number in the dot is even. In this subsection, we present the 
numerical results for the degenerate orbital case. The parameters are chosen as Ac = O.OOSvr, 
Ao = 0.0027r, Aed{= 2t) = 0, U = 0.1 and Jh = for numerical calculation. 

The electron occupation numbers for the even and the odd orbitals increase almost simultaneously 
when the potential ed is dropped as shown in Fig. 7. Small but finite difference of the occupation 
numbers (rig) — (rio) is caused by the difference between Ac and Aq. The two regions, (ric) > (no) 
and (tIc) < (no), are separated at Sd = —0.15 because the effective energy levels of the even and 
the odd orbitals coincide at this point. Conductance at zero temperature, Gp \ is also presented 
in Fig. 7, and is small in all region, especially near the potential ~ —0.15 since (nc) — (no) is 
small. We note that the factor (ne) — (no) in eq. (2.10) is caused by the interference between the 
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Fig. 7. The expectation value of the electron occupation numbers (n^) for the even {p — e) and the odd {p = o) 
orbitals, and the conductance Gp"' for the THi case (n — 1) and the THii case (n = 2) , as a function of the potential 
Ed at T = for the degenerate orbital case. 



tunneling processes via the even and the odd orbitals. The result at T = is contrasted with the 
non-degenerate orbital case shown in § p.3.1| , and is also contrasted with Gp for the case of THii 
which will be discussed in §^. 

Next, we present the temperature dependence of the conductance in Fig. 8. The single particle 
excitation spectra for the potential = —0.15 ((rie) -|- (tIq) = 2.0) are shown in Fig. 9, and 
the Kondo temperature as a function of the potential ea is also shown in Fig. 10. In high 
temperature region (for example T7 = 1.1 x 10^^) of Fig. 8, there are the Coulomb oscillations 
peaks near ~ 0.00, —0.10, —0.20 and —0.30, which correspond to cross points of energies of states 
with different electron numbers; e^i ^ 0, —U, —2U and —3U, respectively. When the temperature 
gradually falls down, the conductance in the potential region —3U ^ ffd ^ once increases, and 
then decreases to approach to Gp ^ at T = 0. This up-and-down behavior of the conductance 
with decreasing temperature seems to reflect the partial growth of the coherency in each orbital 
channel. For example at T = T5, the peak height of the Kondo resonance in even channel is about 
60% of that of the T = limit while 50% in odd channel as seen from Fig. 9. We note the 
relation, = 1.3 x 10"^ ~ = ^K(ed = -0.15) = 1.2 x 10"^. When the temperature decreases 
below the Kondo temperature, the interference cancellation between the even and the odd orbital 
channels becomes gradually complete, and the conductance decreases to approach to Gp''. At 
about T ^ T3 = 1.4 X 10^^ ~ O.IT^, the conductance reaches to the low temperature limit. We 
have weak four peak structures in Gp\ but the peak positions do not coincide with those of the 
Coulomb oscillations. 

We note that the minimum of the Kondo temperature, T^, is much higher than that of the 
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Fig. 8. The temperature dependence of the conductance as a function of the parameter Ed for the THi-2 case. (For 
the degenerate orbital case.) The regions given by broken hnes are the magnetic state and sohd hues region are the 
non- magnetic state (spin singlet state). The suffixes for the temperature are common throughout §2.3.2 and Fig. 
16 in 



non-degenerate orbital case presented in the previous subsection, even though the hybridization 
strengths are same. This behavior has been recognized in the 4-fold degeneracy cases of the impurity 
Anderson model where the parameters satisfy Aq = Aq. 
2.3.3 Effect of Hund's rule coupling 

The energy splitting due to Hund's rule coupling has been observed in the quantum dot which 
have high geometrical symmetry.0 This coupling brings the triplet spin state for the two elec- 
trons configuration of the dot. The Kondo temperature would be strongly reduced because the 
hybridization processes are restricted in this case.© In this subsection the numerical results for the 
presence of Hund's rule coupling are shown. We choose the parameters Ac = O.OOSvr, Aq = 0.0027r, 
Affd = 0, [/ = 0.1 and Ju = 0.020. 

The numerical results of (rip) and Gp^ as a function of the potential are shown in Fig. 11. 
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Fig. 9. The single particle excitation spectra at several temperatures for the potential ea = —0.15 ((rie) = (no) = 
1.0). The results are for the degenerate orbital case. 
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0.05 0.00 -0.05 -0.10 -0.15 -0.20 -0.25 -0.30 -0.35 

Fig. 10. The Kondo temperature as a function of the potential £d for the degenerate orbital case. 



When the potential is dropped, the first electron occupies the dot at ~ 0.00 ((ne) + (no) ~ 
0.5), the second electron at ~ —U + |Jh| = —0.08 {{ric) + (rio) ~ 1.5), the third electron at 
£d ~ -2[/ - I JhI = -0.22 ((nc) + (rio) ~ 2.5), and the forth electron at Ed ~ -3U = -0.30 
{{tLc) + (no) ~ 3.5). (The potential values at Ed = 0, —U + | Jh| , —2U — | Jh| and -3^7 are the 
cross points of the energies of the states with different electron number without the hybridization 
term.) The difference of the occupation numbers (nc) — (no) is brought from the difference of the 
hybridization strengths as previously noted. The way of the electrons occupancy and the behavior 
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Fig. 11. The expectation value of the electron occupation numbers (rip) for the even (p — e) and the odd {p = o) 
orbital, and the conductance Gp"^ for the THi case (n — 1) and the THii case (n = 2), as a function of the potential 
Ed at T = for the presence of Hund's rule coupling case. 



of the conductance as a function of the potential at zero temperature, Gp^ , are roughly similar 



to that of the case shown in § 2.3.2 , in which Hund's rule coupling energy is neglected. But the effect 
of Hund's rule coupling can be recognized in the potential region —2U — | Jh| ^ ed ^ —U + \Jh\ 
(where (ne) + {rio) ~ 2.0). The difference of the occupation number between the even and the odd 



orbitals (rie) — {tIq) is suppressed compared with the results in § 2.3.2 . Therefore the conductance 
at zero temperature, Gp ■* , is very small in the potential region where the interference cancellation 
occurs. 

The temperature dependence of the conductance is presented in Fig. 12. The single particle 
excitation spectra at the potential = —0.15 ((nc) -|- {uq) = 2.0) are shown in Fig. 13, and the 
Kondo temperature as a function of the potential is shown in Fig. 14. 

The Kondo temperature in the region —2U — | Jh| ^ £d ^ —U + | Jh| is far less than that shown 
in Fig. 10. This small Tk is caused by the reduction of the hybridization processes due to Hund's 
rule coupling which restricts the two electrons states in the dot to be the triplet state. 

In high temperature region (for example Tn = 1.1 x 10~^) of Fig. 12, there are the Coulomb 
oscillations at ^ 0.00, —0.08, —0.22 and ~ —0.30. The conductance approaches to Gp at 
about T r^Tsr^ O.IT^ {T^ = TK(ed = -0.15) = 2.1 x 10"^). 

We compare the results in Fig. 12 with those in Fig. 8. The conductance in the two cases 
shows almost identical behavior in the region ^ —U + | Jh| = —0.08 and Sd ^ —2U — | Jh| = 
—0.22. They reach to the lowest temperature behavior for T ^ 10~^ in that region. In the region 
—2U — I Jh| ^ £d ^ —U + \Ju\ in Fig- 12, the conductance initially increases with decreasing 
temperature, and then it decreases in the very low temperature. This behavior is similar to the 
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Fig. 12. The temperature dependence of the conductance as a function of the parameter ed for the THii-3 case. 
(Hund's rule couphng case.) The regions given by broken lines are the magnetic state and solid lines region are the 



non- magnetic state (spin singlet state). The suffixes for the temperature are common throughout §2.3.3 and Fig. 
17 in 



up-and-down behavior in the region of —2U ^ ^ —U in Fig. 8. However in the present case, the 
temperature at which the conductance reaches to the lowest temperature behavior is much lower 
because is very low. The difference of the occupation number is suppressed by Hund's rule 
coupling, so the conductance Gp ^ in this region becomes almost zero. 

The qualitative features of the temperature dependence of the spectra shown in Fig. 13 are 
similar to those shown in Fig. 9 except that the temperature scale of the former is very low. The 
peak at Lij ~ in pe{uj) gradually begins to grow up at about T ~ 10~^, while the peak in /9o(w) 
starts to grow up at about T ~ 10^''. The difference of the temperatures at which the Kondo 
peak starts to grow up in each channel is enhanced by Hund's rule term, though Ac and Aq are 
chosen commonly to those of Fig. 9. But we note that the peaks reach almost simultaneously to 
the low temperature limit at about T ~ 10~^. The finite conductance value in the temperature 
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Fig. 13. The single particle excitation spectra at several temperatures for the potential £d = —0.15 {{uc) + (tiq) — 
2.0). The results are for the presence of Hund's rule coupling case. 
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Fig. 14. The Kondo temperature as a function of the potential Ed for the presence of Hund's rule coupling case. 



region 10~^ ^ T <i 10~^ reflects the growth of the peak in the even channel. The decrease of the 
conductance in 10^^ ^ T ^ 10~^ reflects the growth of the peak in the odd channel. Hund's rule 
coupling prevents to observe the low temperature limit behaviors, but it makes wide an intermediate 
temperature region where a part of the orbital is affected by the Kondo effect. 

§3. Two Conduction Channels in Each Lead 

In we have discussed the situation classifled to the THi that dot's orbitals hybridize only to 
one-dimensional degree of freedom in each lead. In this section we consider also the perpendicular 
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components in the leads as the degree of freedoms. We assume that there are two orbitals in the 
dot, which are even and odd, respectively, under the zx mirror. These are assumed to be even 
under the yz mirror. We name these orbitals in the dot as the even orbital and the odd orbital, 
and denote them by dea and doa- (Note that in §|2|, the even and the odd were denoted with respect 
to the yz plane.) We assume that we have two conduction channels in each lead, which are even 
and odd, respectively, under the zx mirror. These are named as the even and the odd channels 
in the left and the right leads. We write them as ci^ckcn '^hokai ^^,ck<j 

and CRofcfT) respectively. The 

Hamiltonian Hi and -f^i-d is written as follows; 

Hi = ^{ekcl^ka'^Lcka + ekcloka^^Loka) 
ka 

+ Y.(^k4^ekaCRek<7, +efc4oA:aCRofc^), (3-1) 
ka 

Hi-d = -y= '^{Vcdl^ci.cka + KdLcLofcCT + h.c.) 

^ ka 

X!(^edLcRefc(T + K^LcRofca + h.c). (3.2) 

^ ka 

When the operators, Sfco- and a^o- in eqs. (2.5) and (2.6), are replaced respectively by newly defined 
ones, Ska = (cLefcfT + CKeka)/V^ and = (cLofca + CRofca)/\/2, we have the identical form for the 
Hamiltonian. 

But the conductance formula is different from that in §|2[ We have the situation classified to 
the THii that the tunneling via the even and the odd channels are independent in the tunneling 
between leads. The conductance formula is given as follows,© 

= fde(-^).A,p,ie). (3.3) 

p=e,o V / 

Here Pp{£) is the single particle excitation spectrum of the dot orbital with symmetry p. At T = 0, 
the conductance formula is reduced to the following expression. 



gP = ^-^T.-^-'{^M\- (3-4) 



,(2) 

p=e,o 

The numerical calculations are performed for the three cases; (1) non-degenerate orbital case, (2) 
degenerate orbital case, and (3) presence of Hund's rule coupling. The parameters are chosen to be 
same as in §^, for comparison. Then the electronic states of the both systems are identically mapped. 
The conductance for three cases is shown in Figs. 15, 16 and 17, respectively. (Conductance at 
r = has been shown in Figs. 2, 7 and 11, respectively.) As note previously, the interference effect 
of tunneling processes via the even and odd orbitals do not occur in the present THii case. Therefore 
the phenomena due to the interference in THi case do not appear in the present. Maximum of the 
intensity at T = has 2 x 2e^//i, twice of the unitarity limit value of the single conduction channel 
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0.10 0.00 -0.10 -0.20 -0.30 -0.40 

Fig. 15. The temperature dependence of the conductance as a function of the parameter for the THii-1 case. 
(Tunneling processes via the even and the odd orbitals axe independent. For the non-degenerate orbital case.) 
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Fig. 16. The temperature dependence of the conductance as a function of the parameter for the THii-2 case. 
(For the degenerate orbital case.) 

case. In Fig. 15, the differences of the conductance behaviors from those of Fig. 3 are not well 
recognized, because the interference is not so significant in the non-degenerate orbital case. 

In Fig. 16 of the degenerate orbital case, the conductance in high temperature region (for example 
Ty = 1.1 X 10~^) behaves the Coulomb oscillations type similar to that of Fig. 8. As temperature 
decreases, it increases gradually filling the valley between four peaks. At T ~ Ts = 1.3 x 10~^, 
the conductance is rather large indicating the partial growing up of the Kondo resonance in each 
channel. But we do not have the up-and-down behaviors such that shown in Fig. 8. At about 
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Fig. 17. The temperature dependence of the conductance as a function of the parameter Sd for the THii-3 case. 
(Hund's rule couphng case.) 



T ^ T3 = 1.4 X 10"'^ ~ O.IT^, the conductance approaches to the low temperature hmit which has 
a broadened peak structure. In Fig. 17, for Hund's rule coupling case, the conductance behavior 
is similar to that of the degenerate orbital case. But the conductance in the potential region 
£d ~ —0.15 does not increases to the low temperature limit values until the temperature becomes 
very low, T ^ T3 = 2.4 x 10^^ ~ O.IT^. Therefore in the wide temperature range we can see two 
peak structures similar to that of the non-degenerate orbital case. 

§4. Summary and Discussion 

The Kondo effects in the quantum dot systems have been investigated for the cases that (THi) 
there is one conduction channel in each lead and the tunneling processes via the different orbitals 
in the dot interfere each other (in §^), and (THii) there are two conduction channels in each lead 
but the tunneling processes do not interfere (in §^). In both cases, the following three cases for 
the orbital degeneracy have been examined; (1) non-degenerate orbital case, (2) degenerate orbital 
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case, and (3) presence of Hund's rule coupling case. The linear response conductance for the bias 
voltage has been calculated as a function of the gate voltage. The numerical calculations have been 
performed in wide temperature range by using NRG method. 

The Coulomb oscillations are observed in high temperature region commonly in the all cases, but 
various types of conductance behavior appear reflecting the growing up of the Kondo resonance 
at low temperatures. The low temperature limit behaviors appear in the temperature region T ^ 
O.IT^ ~ 0.2T^, where is the minimum value of the Kondo temperature when is changed. 

In the non-degenerate orbital cases, THi-1 and THii-1, the pair of the two peaks of the Coulomb 
oscillations gradually changes one flat peak structure as temperature decreases below 0.2T^. The 
conductance behaviors in these two cases are not so different because the interference effect is not 
significant. 

In the degenerate orbital case of THi-2, the interference cancellation between the tunneling 
processes via the different orbitals causes small conductance at temperatures below O.lTj^. But 
we have the characteristic up-and-down behavior of the conductance with increasing temperature 
in the intermediate temperature range. This is caused by the partial destruction of the coherency 
due to the partial breaking of the Kondo singlet state in each orbital channel. In the case of TIIi-3 
(Hund's rule coupling), the temperature region where the conductance shows the characteristic 
behaviors due to such partial destruction of the coherency shifts to low temperature side, because 
the Kondo temperature is largely reduced. 

In Fig. 16 of TIIii-2 and Fig. 17 of THii-3 in four Coulomb oscillations peaks gradually 
become one broad peak structure as temperature decreases below O.IT^, because the interference 
between two channels does not exist. But we have a regime indicating the partial growing up of 
the Kondo resonance in the intermediate temperature region. In the case of THii-3, we have two 
peak structures similar to that of the non-degenerate orbital case in wide temperature range. 

Electronic states of the cases THi and THii are mapped on with each other, but the tunneling 
current in low temperatures shows quite different behaviors. The sensitivities of the conductance 
to the temperature and to the gate voltage are caused by the difference of the Kondo temperature 
at each potential, are inherent in the interacting systems. 

Here we roughly estimate the characteristic temperature by assuming U ~ lmeV(~ lOK) and 
Ap ~ 0.1meV(~ IK). The Kondo temperature is estimated to be about 3mK (Aq ~ 0.06meV) 
and 30mK (Ag ~ O.lmeV) for the case (1), lOOmK for the case (2), and lO^^mK for the case (3). 
We note, however, the Kondo temperature is very sensitive to the hybridization strength. 

Recently experimental work has been performed with controlling A by Goldhaber-Gordon et 
The change of the conductance due to the Kondo effect seems to appear in their system. The 
results seem to correspond to our numerical results at higher temperature region shown in Fig. 3 
and Fig. 15, 10"^ ^ T ^ 10"^^ qualitatively. 
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Appendix: Non-Degenerate Orbital to Degenerate Orbital 

In and ^, we have used models that the energy separation of the orbitals, Ae^, is extreme. It 
is large (in § 2.3.1 and Fig. 15 in §^) or zero (in § 2.3.2 , 2.3.3 and Fig. 16, 17 in Here we change 
gradually Aej and study when the low temperature properties changes from the non-degenerate 
orbital case to the degenerate orbital case. The parameters for the hybridization strength are 
chosen to satisfy the relations Ac = Aq = A, for simplicity. (For the case of THi, the tunneling 
current between leads disappears by this setting. But in Appendix, we mainly concentrate on the 
effects of the energy separation on the Kondo effect.) The parameters are chosen as A = 0.00257r, 
U = 0.1 and Jh = 0. The potential depth is fixed at = —0.15, therefore the system is in the 
electron-hole symmetry. (The electron occupation number in the dot satisfies {n) = (ric) + (rio) = 2, 
exactly.) 

The electron occupation numbers for the even and the odd orbitals are shown as a function of 
Ae^ in Fig. 18. As Ae^ increases, difference of the occupation number steeply increases at about 




1.0x10 



1 .0x1 0" 



Fig. 18. The electron occupation numbers for the even and the odd orbitals as a function of the orbital energy 
separation Asd- Inset figure is the magnetic excitation spectrum at Ed = 0, T = 0. The Kondo temperature at 
Aed = is estimated to be about T — 1.5 x 10""^ from the peak position of the spectrum. 



Affd ~ 2x 10 ^ (where (ne) ~ 1.5 and (rio) ~ 0.5). The Kondo temperature at Ae^ = is estimated 
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to be about Tk = 1.5 x 10^^ from the magnetic excitation spectrum shown in the inset of Fig. 18. 
We may conclude that the cross over occurs at about Ae^ ~ TK(Aed = 0). We note that Ae^ is 
much smaller than the hybridization width, A. 
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